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Abstract 

We study the transverse single spin asymmetry in rye, 6 production in polarized hadron collisions, 
employing the collinear twist-3 approach in combination with the color singlet model. Our main 
focus lies on the contribution from the twist-3 Efremov-Teryaev-Qiu-Sterman function. By ex- 
trapolating the derived spin dependent cross section to the small transverse momentum region, 
consistency between the collinear twist-3 approach and the transverse momentum dependent fac- 
torization approach is confirmed. As a byproduct of this work, we identify a term contributing to 
the scale evolution of tri-gluon correlations in the flavor-singlet case which was previously missed. 
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I. INTRODUCTION 



Attempts to understand the large transverse single spin asymmetries (SSAs) observed in 
the various high energy scattering processe presently form an active field of research. The 
experimental observations arc a great theoretical challenge as the coUinear leading twist 
contribution to the asymmetries is proven to be proportional to the quark mass [1, 2] and 
thus very small. However, during the past few decades, remarkable theoretical progresses 
was achieved by following mainly two approaches: one approach is based on transverse 
momentum dependent (TMD) factorization [3-8] and the other on coUinear twist-3 factor- 
ization [9-17]. In TMD factorization, a naive time reversal odd TMD distribution, known as 
Sivers function [3] describing the correlation between the parton intrinsic transverse momen- 
tum and hadron transverse spin vector is responsible for the asymmetries. In the coUinear 
twist-3 approach the SSAs arise from twist-3 quark gluon correlators so-called Efremov- 
Teryaev-Qiu-Sterman functions (ETQS) Tp{x,x) [9, 10]. These two framework have their 
own kinematic regions of validity and were shown to produce the same results in the overlap 
of these kinematic regions [14-16, 18]. 

The fundamental property of both mechanisms is that the imaginary part necessary for 
non- vanishing SSAs is dynamically generated through initial/finial state interactions. In the 
case of TMD factorization, the effect of these initial/finial state interactions is encoded in the 
process dependent gauge link appearing the matrix element defining the TMD distributions. 
Due to different structures of these gauge links in different processes, naive time reversal 
odd TMD distributions, like the Sivers function possess a very unique, modified universality 
property, namely, that the Sivers functions in DIS and Drell-Yan process differ by a minus 
sign [6, 7]. 

While SSAs in DIS have been measured by HERMES [19] and COMPASS [20], a polarized 
Drell-Yan measurement is not yet available, such that this prediction could not yet be tested. 
Alternatively, one can simply fit the ETQS function related to the kr moment of the Sivers 
function [22] using input from HERMES and COMPASS [21] and compare it with a direct 
extraction of Tp from the SSA obsercvd for pion production in polarized pp collisions [13, 23- 
25]. Promising early results along these lines, have, unfortunately been caused by a sign 
error, such that the obserevd discrepancies are now referred to as "sign mismatch" of these 
two processes [26]. Our recent work [27] also indicated that this procedure gives the wrong 
sign for SSA in inclusive DIS off a neutron target (for a complete theoretical treatment of 
this asymmetry, see Ref. [28]) if the quark gluon correlation Tp fitted to data for the process 
P t P ~^ ttX were used as the input for the twist-3 quark photon correlation. This suggests 
that the SSAs observed in process like p t P ~^ iiot mainly caused by the Sivers 

mechanism as described in the twist-3 coUinear approach [11, 13]. A possible explanation is 
suggested by the observation [29] that the Collins effect [5, 30] describing SSA as caused by 
parton fragmentation can contribute significantly. This observation has been confirmed by 
a more recent analysis [31]. 

In order to better understand the Sivers mechanism in hadronic collisions, it is desirable 
to investigate SSA for cases of particle production in polarized pp collisions for which the 
Collins effect is absent. Possible options are the SSA in direct photon or jet production 
production p ^ p ^ 'y{Jet)X [10, 32]. Another option which was initially proposed in 
Ref. [33] is the SSA in heavy quarkonium production. The author of the paper [33] analyzed 
SSA in heavy quarkonium production following the general arguments of non-relativistic 
QCD(NRQCD) [34, 35]. In the framework of NRQCD, the heavy quarkonium can be pro- 
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duced at short distance not only in the color singlet but also the color octet configuration. 
The transition from heavy quark and antiquark pair to a quarkonium state is treated as 
a non-perturbative process and encoded in the long distance universal matrix elements 
which are characterized according to the velocity expansion of NRQCD. As demonstrated 
in Ref. [33], the SSA in heavy quarkonium production offers a unique way to investigate its 
production mechanisms since the SSA crucially depends on the final state interactions which 
differ significantly for heavy quarkonium production in the color singlet and the color octet 
channel. Recently, a SSA in p t p ^ J/ ipX has been measured at PHENIX and found to 
be sizeable [36] . According to the analysis of Ref. [33] , this hints at a dominance of a color 
singlet mechanism at low transverse momentum and at a non-zero gluon Sivers effect. Also, 
in Ref. [37], the color evaporation model in combination with TMD factorization was used 
to study the SSA for J/tp production in ep t collisions. 

Following this research line, we apply the coUinear twist-3 approach to compute the 
SSA for heavy quarkonium production in polarized pp collisions. In the present work, the 
earlier pioneering analysis [33] is extended and refined in the sense that we compute the 
transverse momentum dependent behavior of the spin asymmetry, and also take into account 
the hard gluon pole contribution. We only focus on r]c,h production in this paper, though our 
formalism can be easily extended to SSAs for the production of other C-even quarkonia, such 

XcOi XfeOi Xc2 and Xb2- In leading order, a C-even quarkonium can be produced through 
two gluon fusion, which makes it a promising way to access both unpolarizcd and linearly 
polarized gluon TMD distributions [38-42] . The application of TMD factorization in these 
processes has been justified by a recent NLO calculation [43]. Experiments measuring such 
flc,b asymmetries could be performed at RHIC and a proposed fixed target experiment at 
LHC (AFTER) [38-40]. LHC would allow us to study the gluon polarization effect in a deep 
saturation regime[44] due to the very high energy of LHC and relative low mass of C-even 
charmonium. In Refs. [45, 46], the promising channels to detect r]h were discussed. 

The analysis [34] based on NRQCD suggested that for C-even quarkonium production, 
the color octet contribution is suppressed while the color singlet contribution dominates, 
in particular at low transverse momentum [47], though complications may arise at large 
transverse momentum [48, 49]. Moreover, C-even bottomonium production in the color 
octet configuration can certainly be neglected [38, 39, 45], since it is strongly suppressed 
in the velocity expansion. On the other hand, the SSA in C-even heavy quarkonium pro- 
duction receives a contribution not only from ETQS function but also from the tri-gluon 
correlation [50]. Nevertheless, we restrict ourselves to color singlet production, i.e. the Tp 
contribution, and leave the color octet part and tri-gluon correlation for a future study, as we 
are primarily aiming at establishing a formalism that combines coUinear twist-3 techniques 
and NRQCD factorization to describe SSA in heavy quarkonium production in the present 
paper. 

Our calculation is carried out in the covariant gauge closely following the techniques out- 
lined in Refs. [10, 11, 13] in which the SSAs for pion and direct photon production in hadronic 
collisions were computed. Recently, in analogy to the SSAs, the double spin asymmetry Alt 
in the same processes was also studied within the coUinear twist-3 framework [51]. Unlike 
the SSA in pion or photon production, TMD factorization could be applied in the kine- 
matical region where the transverse momentum of the heavy quarkonium is much smaller 
than its mass. We extrapolate the complete coUinear twist-3 result to the small transverse 
momentum region and find that it is consistent with the result obtained from TMD fac- 
torization. The key step in establishing this connection is to derive the large transverse 
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FIG. 1: Diagrams contributing to the unpolarized cross section at order O(a^) and 0{a^). Gluon- 
gluon scattering at order O(q^) is not shown. Diagrams with permutations of the gluon hnes are 
also not shown. 

momentum gluon Sivers function and relate it to the ETQS function. Doing so, we find that 
the hard gluon pole contribution to the gluon Sivers function was overlooked in the previous 
literature [33]. We will also comment on the scale evolution of tri-gluon correlations for the 
flavor singlet case in the following section. 

The paper is organized as follows: in the section 2, we first briefly review ?7c,b production 
in unpolarized hadron collisions. Next we derive the spin dependent differential cross section 
for rjc^b production in hadronic collisions, presenting the result in the limits of high and low 
transverse momentum. In particular, we discuss the matching between TMD factorization 
and the coUinear twist-3 approach at small transverse momentum. The paper is summarized 
in section 3. 

II. CALCULATION OF THE UNPOLARIZED AND POLARIZED CROSS SEC- 
TIONS 

We start by introducing the relevant kinematical variables. For the process under con- 
sideration. 



we define the 4-momenta and polarization vector of the incoming nucleons A, B and outgoing 
heavy quarkonium rjc^b as indicated. The Mandelstam variables are S* = (P + P')^, T = 
[P' — lY and U = {P — ly . The corresponding Mandelstam variables on the partonic level 
are given by s = {xP+x' P'Y , t = [x'P'—lY and u = {xP—l)"^ where x, x' are the longitudinal 
momentum fractions carried by the partons from nucleon A and B respectively. The squared 
invariant mass of the heavy quarkonium is = and M = 2Mq up to small relativistic 
corrections, where Mq is the heavy quark mass. In the unpolarized case, to lowest order 
0(q;^), one only has the gluon fusion process shown in Fig. 1(a). This production mechanism 
dominates at small transverse momentum where the mean produced quarkonium transverse 
momentum is of order Aqcd- At order O(a^), there are much more contributions from the 
two- by-two scattering processes illustrated in Fig. 1(b) and Fig. 1(c). In these processes, 
the large heavy quarkonium transverse momentum is generated by recoiling against an 
outgoing gluon or quark. Note that heavy quarkonium production at O(a^) from gluon 
gluon scattering is not shown in the Fig.l. 




(1) 
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Within the so-called color singlet model [52-54] , the amphtude for the processes discussed 
has the general form, 



M 



J (0iTr[O(U)0(/,g)] (2) 



where 2q is the relative momentum of the heavy quarks. Here, q) is the Bethe-Salpeter 
wave function of the produced bound state. 0(1, q) is the perturbativc part of the diagrams 
shown in Fig.l with the heavy quark legs cut off. Due to the non-relativistic nature of heavy 
quarkonium, the relative momentum is assumed to be small with respect to the quark mass 
M/2. For the production of S-wave states, one can simply neglect the dependence of 0{l, q) 
on the relative momentum. This leads to, 

^R,momvss.m (3) 

V47r 



Here, -Ro(O) is the value of the radial S-wave function at the origin, and VsSz (0) is the spin 
projection operator. According to the color singlet model, the heavy quark pair produced in 
the partonic scattering can evolve into a quarkonium non-pcrturbatively only if they have 
the same quantum numbers as the corresponding quarkonium. This implies that the spin 
projector operator associated with r]c,h production has the form. 

With these calculational recipes, one derives the unpolarized differential cross section for 
the gluon fusion channel [54], 

-aoG,{z)G,{z') (5) 



dydH 



where y is the rapidity of the produced bound state. Gi{z),Gi{z') are the unpolarized 
collinear gluon distribution functions of two incoming hadrons. The longitudinal momentum 
fractions are constrained by kinematics: z = ^e^, z' = In Eq.(5), ctq is given by. 



^2 r2 „ ,2 



For the quark gluon scattering channel, we have [54] , 

where fi{x) is the unpolarized collinear quark distribution. The index a runs over all quark 
flavors. For the qq channel, the unpolarized differential cross section reads [54] , 

da -8 as ,_^2V- f dx dx' - {s-M^f-2iu 



dyd'^lj 



^ — 8 tts , ^9 f dx dx . ,. (s — M ) — 2iu ^, ^ - ^ , ^9x / x 



One notices that apart from the color factor, the hard coefficient appearing in the qq channel 
can be obtained from those of the quark gluon channel by crossing £ s. The result for 
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FIG. 2: Sample diagrams giving rise to the SSA in hadronic r/c;, production from quark gluon 
channel. 

the gluon-gluon scattering channel can also be found in Ref. [54]. Next to leading order 
corrections to C-even hadronic quarkonium production have been calculated in Refs.[55, 56]. 
Now let's proceed to derive the polarized cross section. To generate the spin asymmetry, 
one additional gluon must be exchanged between the active partons and the remanent part 
of the nucleon target as shown in Fig. 2 and Fig. 3. The hard part, if an additional gluon 
is attached, can be calculated perturbatively, while the non-perturbative part describes 
the relevant three parton correlations. As stated in the introduction, we only focus on 
the contribution in which the quark-gluon correlation in the transversely polarized nucleon 
enters. This is precisely the ETQS function which is defined as. 



where we have suppressed Wilson lines and have indicated the nucleon mass by . We also 
introduced the lightcone vector n = (l"*", 0~, 0_l), whose conjugate vector is n = (0+, 1~, 0_l). 
Note that our definition of the ETQS function differs by a factor 27rMjv from the conventions 
used in Refs. [14, 15]. This ETQS function plays an important role for SSA phenomenology. 
Its scale evolution has been derived in Refs. [57-64]. 

Similar to the SSA in the Drell-Yan process, the strong interaction phase factor neces- 
sary for having a non-vanishing spin asymmetry arises from the interference between an 
imaginary part of the partonic scattering amplitude with an extra gluon, as shown in Fig. 2 
for the quark-gluon scattering channel, and the real scattering amplitude without a gluon 
attachment. The imaginary part is due to the pole of the parton propagator associated 
with the integration over the gluon momentum fraction Xg. This effectively implies that 
one of the internal parton lines goes on shell. To isolate the imaginary part of such poles, 
the distribution identity: = PV^ ^in6{x) was used. Depending on which propagator's 
pole contributes, the amplitude may get contributions from Xg = { "soft-pole" ) and Xg 
("hard-pole" ) [14, 15]. Both types of gluon poles show up in our calculation. 

As the reader can find the relevant technical details for our calculation in the litera- 
ture [10-15] we sketch here only some key steps of such twist-3 calculations. 

As mentioned before, we carry out the calculation in the covariant gauge, in which the 
leading contribution of the exchanged gluon is the "plus" component A'^. The gluon's 
momentum is dominated by the XgP + kg±, where Xg is the longitudinal momentum fraction 




(9) 



6 



(6) 



FIG. 3: Sample diagrams giving rise to the SSA in hadronic r/cfe production from the quark- 
antiquark channel. 



with respect to the polarized proton. In order to calculate consistently with twist-3 accuracy, 
one has to expand the hard parts in the gluon transverse momentum and keep the terms 
linear in kgj_. Then the kgj_ factor can be combined with to yield d-^A~^, which is an 
element of the field strength tensor F^^. After adding the term proportional to d~^A-^ of the 
same tensor [12, 65], the soft part can be rewritten in the form of the ETQS function. Making 
use of the ingredients described above, the calculation is straightforward. However, the 
derived complete results for the spin dependent cross sections is rather length. To highlight 
some interesting features of our results, in the following, we will only present the expression 
in the limit of large transverse momentum (/^ 3> M) and small transverse momentum 
(Aqcd <^ ^± ^ M) instead of the full expression. At large transverse momentum, after 
dropping all terms suppressed by powers of M/l±, one ends up with the polarized differential 
cross section for the qg channel. 



da 



dyd'^lj 



X 



2TUx,x) - X 



d_ 

dx 



Tj^ [x . X , 



Tp(^x, x) 

and for the qq channel 
da 



2s^ + is 1 



t3 



Nc2ut 



X x' —i 

~2 '^f 

Tp \x^x 



u) 



s - t 



Nc + sH - t 
~2 ^ 



f'3 



(10) 



dyd'^lj 



X 



8\ ^ , etc , ^9 , ^ f dx dx' 1 77^,,,,,^ 



u 



2T^{x,x)-x[-^T^{x,x] 



-1 e + e 

2Nr -s^ 



Tp(^x, x) 



iVc 2P + is 1 i 



Nc2us 



Tp I X, X 



-t \ Nc 2t^ + t'^s - s- 



s-t 



—utp 



A few remarks on these analytical results are in order. First, as usual, the SSA is suppressed 
by the factor Mm/1± at large transverse momentum. Second, one notices that the derivative 
term and non- derivative term from the soft gluon pole contribution can not be combined 
into a compact form as it can be done for the polarized cross section for pion production 
in hadron collisions [13]. Also, let us note that the hard gluon pole contribution survives 
even though the heavy quarkonium mass is neglected, in contrast to the cases of pion and 
direct photon production for which the hard gluon pole contribution is absent. Finally, it is 
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(a) 



FIG. 4: Diagrams contributing to the gluon Sivers function at large transverse momentum in the 
flavor singlet case. Fig. 4(a): soft gluon pole contribution; Fig. 4(b): hard gluon pole contribution. 



observed that except for the different color factors, the soft gluon contribution in the quark 
gluon channel can be obtained from the quark antiquark channel by crossing s ^ t, while 
the hard coefficients associated with the hard gluon pole contributions for the qg and qq 
channels differ by a minus sign after crossing s -H- 1. 

Now let us discuss our results in the low transverse moment limit. In this limit, we 
neglect all terms suppressed by the power of l']_/M'^ while keeping those enhanced by the 
factor In order to extrapolate our results to the small transverse momentum region, 

one has first to expand the delta function in l±, 



1 



(12) 



where ^ = z/x and ^' = z' /x'. The 'plus' prescriptions is defined in the standard way. With 
this expansion, the spin dependent cross section in the kinematical region Aqcd ^ l± ^ M 
is given by. 



da _ Nc as Mn ^ , ^ , n f dx 



S±^.huGi{z') j 



dyd'^l± 2 n l\_ J X 

X I (x-^Tpix, x)] [1 + (1 - 01 - Tpix, x)2{l - 0^ - Tf{x, x - z'^ ~ ^ 



=^0^:^ 7^^'^ S^^l^uGi{z) / — <Tf{x,x) Tf{x,x - z)——\ {13) 

2 -K i\ J X y i i ] 

In the second step of the derivation of above formula, we carried out the integration over 
X by parts. It is interesting to note that the hard coefficient associated with the soft gluon 
matrix element is the well known splitting kernel Vgq. Here, it is worthwhile to point out 
that no leading power contribution at small l±_ comes from the quark-antiquark channel. 

On the other hand, when l±_ <^ M, transverse momentum dependent factorization can 
be applied. In the TMD factorization approach, the SSA in 77^6 production is generated by 
the gluon Sivers function Gj-j^ through the gluon-gluon fusion channel, 

^o ^'^^y^" I d'k^d'k'^^^6\h -k^- k'^)G,{z',k'^)Gi^^y{z,k^) (14) 
Mn j I, 



dydHj 
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where Gj-j^ ^yi^y ^'^) denotes the gluon Sivers function. The subscript "DY" indicates that 
the gluon Sivers function contains a past-pointing gauge hnk built up through initial state 
interactions, similar to that in the Drell-Yan process. 

When k± is of the order of Aqcd, the gluon Sivers function is an entirely non-perturbative 
object. However, in the kinematic region Aqqd -C <C M, TMD factorization still holds 
and at the same time the function G^'-^^^y (2;, A;^) can be calculated in terms of twist-three 
parton correlation function within perturbative QCD. It can receive contributions from both 
the ETQS function and tri-gluon correlations. In the current case, we are interested only 
in the former one. Our perturbative calculation follows a similar procedure as in [14, 15] 
resulting in 

GlW^.y = f /f {.(...)i±(lzi)!_T,(...-.)^} (15) 

where the soft gluon pole contribution is generated by the diagram shown in the Fig.4(a), 
and the hard gluon pole contribution arises from Fig.4(b). This result was first derived 
in Ref. [33], though the second term was missed there. At this point, we would like to 
mention that the hard gluon pole term also contributes to the scale evolution of the tri-gluon 
correlation for the flavor singlet case, which was overlooked in the literature previously [57, 
60] ^. In addition, due to the existence of diagram Fig.4(b), the chiral partner of the ETQS 
function, Tf{x,xi) (notation used in Ref. [16]) also contributes to the scale evolution of 
tri-gluon correlations [67]. 

By inserting Eq.l5 into Eq.l4 and making the approximation — k± — k'j_) ^ S'^{1± — 
k±), we reproduce the result obtained in the coUinear twist-3 approach. Therefore, for the 
observable under consideration, we have obtained a unified picture in the kinematical region 
where TMD factorization and the coUinear twist-3 approach both apply, as was found in 
many other cases [14-16]. 

To end, let us briefly discuss a possible future numerical study. In order to estimate 
the size of SSAs for ric,b production, we must determine the input for the ETQS functions 
including the diagonal and off-diagonal contributions, as the polarized cross section depends 
on the both, soft and hard gluon pole contributions. Unfortunately, the off-diagonal con- 
tributions to Tp needed for this SSA observable are not as well-determined as the diagonal 
pieces which can be related to moments of the quark Sivers function. In Ref. [57], a Gaussian 
form was assumed for Tp{x,Xi) with a maximum at x — Xi. This study was done in the 
context of the evolution of Tp{x,x). In Ref. [66], an analysis of higher- twist functions was 
conducted using light-cone wave functions that include qqqg Fock states. In contrast to [57], 
this study found that Tf{x,Xi) reaches its maximum when x ^ Xi and some of its lowest 
values when x = xi. These completely different behaviours clearly demonstrate our lack of 
knowledge for Tf{x, xi). Nevertheless, it would be useful to determine the potential impact 
of Tf{x, Xi) on the size of SSAs in rjcb hadronic production. 



This term is present in [60] in the evolution equation for arbitrary gluon momenta, but missing in their 
answer for the gluon pole limit. 
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III. SUMMARY 



In summary, we have calculated the transverse single spin asymmetry in hadronic 7]c,b 
production by employing the coUinear twist-3 approach in combination with the color sin- 
glet model. We discussed the behavior of the spin asymmetry at high and low transverse 
momentum. In particular, at low transverse momentum, a match between the TMD fac- 
torization approach and the coUinear twist-3 formalism has been found after re-examining 
the derivation of the gluon Sivers function at large transverse momentum. As a byprod- 
uct of this work, we identified a term, which was missed so far, that contributes also to 
the scale evolution of the tri-gluon correlation for the flavor singlet case. In addition, we 
have briefly outlined a plan for a future numerical studies. Let us mention that it would 
be feasible to measure this asymmetry at RHIC and a proposed flxed target experiment at 
LHC(AFTER) [38-40]. 

We emphasize that the transverse spin physics and the heavy quarkonium production 
physics could mutually beneflt from studying this observable as explained in the following. 
First of all, due to the sign mismatch issue [26] , a doubt on the validity of the coUinear twist- 
3 approach in hadron collisions has arisen. SSA in heavy quarkonium production provides a 
clean way to test this approach, as compared to that in pion production where the Collins 
mechanism could dominate the asymmetry [29]. Furthermore, one notices that the spin 
asymmetry is independent of the radial wave function at the origin Rq and only sensitive to 
the heavy quarkonium production mechanism. Therefore, measuring this observable would 
provide us an unique chance to pin down the ratio between the values of the color singlet 
long distance matrix clement and the color octet matrix element. 

There are a number of directions in which our work could be extended. First, the color 
octet contribution to SSA should be taken into account as it plays a role at large transverse 
momentum. Second, it is natural and straightforward to study the SSAs for other C-even 
heavy quarkonium production in the framework outlined in this paper. It is also possible 
to study SSAs for hadronic J/ip^T production. To do so, we have to compute the SSAs 
generated from the tri-gluon correlation since it is the only contribution appears at the 
non-trivial leading order in J/ip^T production. Finally, one can calculate the SSAs for 
photonic/electronic J/'0,T production using the coUinear twist-3 approach in combination 
with the NRQCD framework. These would be the relevant observables at the future EIC [68, 
69]. 
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